Superconductivity in bundles of a mixture of doped carbon nanotubes 
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Using inhomogeneous Bogoliubov-de Gennes formalism we study superconducting properties of 
bundles of single wall carbon nanotubes, consisting of a mixture of metallic and semiconducting 
nanotubes, having different critical transition temperatures. We investigate how the averaged su- 
perconducting order parameter and the critical transition temperature depend on the fraction of 
semiconducting carbon nanotubes in the bundle. 



Single wall carbon nanotubes (SWCNT) represent a 
unique class of quasi-one dimensional nanoscale systems, 
exhibiting various interesting phenomena. Among other 
exciting features, it was demonstrated that individual 
single wall carbon nanotubes may have intrinsic super- 
conducting properties^. However, because of their ex- 
tremely small diameter (just few nanometers), and thus 
strongly one dimensional character, the superconduct- 
ing order parameter may have significant "phase slips" 
due to thermal and quantum fluctuations, leading to 
a finite conductivity in the system below the critical 
temperature^i^. Carbon nanotubes can form bundles and 
ropes^, with tens and hundreds of individual SWCNTs 
in the bundle, coupled to each other by dispersive Van 
der Waals forces. Such kind of system may exhibit re- 
duced "phase slips" effects, due to three-dimensional cou- 
pling of the nanotubes in the bundle and as a result, 
much stronger conductivity drop below the critical tem- 
perature. The overall length of a SWCNT in the bun- 
dle also plays a significant role. For example, reducing 
the bundle's length to 300 nm destroys the supercon- 
ductivity in the system due to increasingly high quan- 
tum fluctuations^. Generally speaking, for nanoscale sys- 
tems with the quantum level spacing approaching the 
superconducting gap energy A, the superconductivity 
vanishes^. 

It is expected that doping of SWCNTs in a bundle by, 
for example, boron, may significantly improve their su- 
perconducting properties. It is believed that a proper 
level of doping may result in the Fermi level at a one 
dimensional singularity of the energy spectrum and thus 
in a higher density of states (DOS), that will lead to a 
higher critical temperature T^. In particular, we suggest 
here that such kind of mechanism of doping enhanced Tc 
may be much better pronounced in the case of semicon- 
ducting SWCNTs, which may have higher DOS due to 
lower in energy van Hove singularities. This is in con- 
trast to metallic SWCNTs, where singularities in the 
DOS are much higher in energy, and start being filled 
much later during the doping process (according to the 
famous Kataura plot^). Therefore, a bundle consisting of 
doped semiconducting nanotubes could be a much better 
superconductor, compared to a bundle made of metallic 



SWCNTs. 

However, synthesis of SWCNTs by all currently known 
methods results in a mixture of semiconducting and 
metallic nanotubes. Since the nanotubes after the syn- 
thesis initially are not doped (or unintentionally slightly 
p-type doped, e.g. by oxygen of atmosphere), those are 
only metallic tubes, which may have superconducting 
transition, while semiconducting tubes will be "diluting" 
superconductivity in the bundle by the inverse proximity 
effect^. Upon doping (i.e. by electrochemical charging), 
the semiconducting tubes can become superconducting 
with a higher superconducting gap and thus a higher 
than in metallic nanotubes. 

Therefore, one should be able to estimate spatially av- 
eraged order parameter and the corresponding effective 
critical temperature for a bundle consisting of a mixture 
of SWCNTs of these two model types. From an exper- 
imentalist's point of view it is even more important to 
solve the inverse problem: for a given fraction a of semi- 
conductor SWCNT in the bundle and the experimentally 
determined critical temperature Tc(a), to estimate the 
critical temperature for a bundle, consisting only of semi- 
conductor SWCNTs Tc{a — 1)? It will be also interesting 
to know, can one obtain much higher than in other car- 
bon based nanostructures, and particularly higher than 
in alkali metal doped fullernes. 

Spatial variations of the superconducting order pa- 
rameter are significant for nanoscale systems, including 
nanotubes®'^. In this work we use a microscopic theory 
based on inhomogeneous Bogoliubov-de Gennes equa- 
tions to establish how the superconducting properties of 
a bundle depend on the fraction of doped semiconductor 
nanotubes, with a higher SC order parameter. We as- 
sume that the nanotubes in the bundle are approximately 
of the same radii and tightly packed making a triangular 
lattice in the bundle's transverse section, with the primi- 
tive vectors oi = Rx, 02 = Rx /2+\/^/2Ry. Here af, y are 
the unit basis vectors, and R is the intertube distance. 
The lattice can enumerated by indexes (i,j), which cor- 
respond to the position of a nanotube — dii + a2j, 
but in this work we prefer to enumerate nanotubes in a 
N X N bundle using a single index through the mapping 
k — i+jN, where TV is the number of nanotubes in a raw. 
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In the bundle semiconductor nanotubes are assumed to 
occupy the fraction a of the sites, and metalhc nanotubes 
1 — a, accordingly. In our model the conduction electrons 
can freely travel along the nanotubes, and in this picture 
it corresponds for an electron staying at the same lattice 
site. However, electrons can also hop to the neighbor- 
ing nanotubes (sites). In principle, there may be three 
different hoping constants, with the hopping matrix ele- 
ments t'^'^ equal to either tmm,tms,tss, corresponding to 
the hopping between metallic-metallic (mm), metallic- 
semiconducting (ms), or semiconducting-semiconducting 
(ss) nanotubes. Moreover, these parameters may signif- 
icantly fluctuate from one site to another, due to mis- 
match between SWCNTs of different chirality. In the 
superconducting regime Cooper pairs can be formed and 
can freely move along the nanotubes, and can also hope 
from one tube to another. 

For the description of the system we utilized a tight- 
binding Hubbard Hamiltonian of the form: 

Y ^^m^(l■^)"4.(l■^)'^t(l■^) 

<ri> 

+ E Ct".(r.)".'(r,), (1) 

where a quantum- mechanical operator cj. „ creates an 
electron on site i (using single indexing), the operator 
Crj,CT eliminates an electron from the site j, and ncy{vi) ~ 
cJ._^Crj,o- represents the electron density on site i with 
the spin polarization a. The electron spin, <t, can point 
up or down. Ul^^^ is on site interaction potential. This 
term in a case of attractive interaction Ul^^ < may 
lead to pairing in the nanotube i. V-^^ is a strength of 
the coupling between electrons localized at neighboring 
tubes i and j. 

Using the Boguliubov transformation, which diago- 
nalizes the Hamiltonian Eq.([T]), we arrive to inhomoge- 
neous Bogoliubov-de Gennes equations for the quasipar- 
ticle amplitudes on the lattice i sites ('u„(ri), w„(ri))iS: 




where the kinetic operator ^ and the superconducting 
order parameter A can be represented as: 

iun{ri) = - E^*'^Wn(ri + + (y'(ri) - /Li)u„(ri), 
s 

kvn{Yi) = E ^s{'rt)VniYi ^8) + (rj)i;„ (r^ ) , (3) 
s 

where 8 are the nearest neighbor vectors for a triangular 
lattice, y (r^) is the mean- field (Hartree) potential, is 
the chemical potential. Ag is the conventional, s-type 
order parameter. One should solve Eq.® together with 



the self-consistency conditions: 

A,(rO=E^K(^»+^K(^^) + 

n 

u„(r,)<(r, + 8)) tanh(£;„/2fcBT), (4) 

where the pairing strength V^^^j may depend on the type 
of CN at sites i and j. The s-type order parameter 
(within a given nanotube i) is simply 

A.(r.)=E%^K(rO<(r.) + 

n 

u„(r,X(r,)) tanh(S„/2fcsT), (5) 

Note, the summation in Eqs. ()4l5|) is over the positive 
eigenvalues En only. 

Here we adopted a simplified picture assuming the 
same constant hopping parameter t between any type 
of nanotubes. In this work we considered for simplic- 
ity that the pairing may happen between electrons in the 
same nanotube, therefore neglecting much weaker pairing 
mechanism between neighboring nanotubes. In principle, 
a weak attraction mechanism may stimulate the forma- 
tion of a Cooper pair with one electron in one nanotube, 
and the second electron in one of its nearest neighbors. 
This may result in co- existing order parameters in the 
system, one of the order parameters with the conven- 
tional (s-type) symmetry, and another with unusual sym- 
metry. The co-existence of order parameters with differ- 
ent symmetries, was studied, for example, for Uranium- 
based superconducting materials^-. The possibility of 
such pairing on a triangular lattice may result in uncon- 
ventional superconducting properties. For example, a 2D 
triangular lattice was recently considered as a test-bed for 
a possibility of f-wave spin-triplet superconductivitytl. 

The amplitudes u„(ri), w„(ri) obey the constraints 
/ (ir(|M„(ri)p -f |?j„(ri)p) = 1 for any n (normalization) 
^'^d i]„(|Mn(rOP -f |w„(ri)p) = 1 for any i, i being the 
site index of the triangular lattice. 
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Fraction of semiconductor CN 

FIG. 1. Spacially averaged superconducting order parameter 
< A > (in units of t) as a function of fraction of semiconductor 
nanotubes a for different temperatures. 

We studied how a spatially averaged superconducting 
order parameter < A >= As(rfc) depends on the 
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FIG. 2. Spacially averaged superconducting order parameter 
< A > (in units of t) as a function of temperature for different 
values of a (fraction of semiconductor CN). 
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Fraction of semiconducting CN 

FIG. 3. Critical temperature Tc as a function of the fraction of 
semiconductor nanotubes a. Note a steeper slope for a < 0.5 
(the percolation limit in two dimensions on the triangular 
lattice) . 



fraction of semiconductor SWCNTs a in a iV x bundle 
at different temperatures. For this purpose we gener- 
ated P = 50 realizations for a given number of randomly 
placed semiconductor nanotubes in the bundle. The rest 
of nanotubes in the bundle are assumed to be metallic. 
To model different superconducting pairing strength for 
different types of nanotubes we set = 2t for semi- 
conducting nanotubes and ?7*„j = 0.68i for metallic ones, 
and assumed /i = (half filled band). In our simulations 
we considered 16 x 16 nunotubes in the bundle, forming 
a triangular lattice. 

The results of calculations are shown in Fig.l. At T ~ 
the order parameter scales approximately as a square 
root of the fraction of the semiconducting SWCNTs in 



FIG. 4. Spatial distribution of the superconducting order 
parameter (in units of t) at the percolation regime (a — 0.5) 
and zero temperature. Red dots mark the triangular lattice. 



the bundle cx {T^'-T^)y^+T^. Note the convexity of 
the dependence. However, at finite temperatures, which 
are between the critical temperature of a pure metallic 
and pure semiconducting bundles, the averaged 
order parameter vanishes much faster with the decreasing 
of a. For example, at T = 0.32t, which is close to the 
critical temperature of a pure semiconducting SWCNT 
bundle « 0.35t, the order parameter decreases in the 
exponential fashion with the decreasing of a, and almost 
vanishes at a « 0.5. 

In Fig. 2 we also plot < A > as a function of tempera- 
ture for several values of a. One can clearly see how the 
order parameter vanishes above the critical temperature. 
Note that with the lowering of the fraction of semicon- 
ducting nanotubes a, the temperature dependence of the 
order parameter shows less pronounced phase transition 
because of the "dirty" nature of inhomogeneous spatial 
distribution of the pairing properties, similar to dirty su- 
perconducting transition in case of large concentration of 
impurities and in alloys. 

Using the data plotted in Fig. 2, we calculated how 
the critical temperature Tc depends on the concentration 
of semiconducting SWCNT. We used t = 4.8 x IQ-^eV 
to fit the data in [22-23], so a = 0.6 will correspond to 
Tc ~ 15K. In Figure 3 one can see that the critical tem- 
perature decreases for a < 0.5 with a steeper slope, where 
a = 0.5 corresponds to the percolation threshold in 2D 
triangular lattices (please see Fig. 4). For a < 0.5 weaker 
superconducting metallic CN are arranged in bigger size 
islands. If one would take into account the phase fluc- 
tuations of the order parameter, this decrease of critical 
temperature would be even steeper, because of enhanced 
phase slips in relatively well isolated semiconducting nan- 
otubes. 

It should be noted that our model has general appli- 
cability to any system in which there are two types of 
nanotubes, (or very thin nanowires) with different super- 
conducting pairing strength are coupled in bundles. So 
it also describes the most common case of undoped pris- 
tine SWCNT bundles, which contain 30-40 % of metal- 
lic SWCNT and the rest are non-doped semiconducting 
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SWCNTs, which usuahy do not have carriers. As has 
been shown this bundles have typical of 0.55 K [20- 
21], which according to our model is suppressed by the 
inverse proximity effect from non-superconducting un- 
doped semiconducting tubes. According to our model if 
100 % separated only metallic tubes are in the bundles, 
then the gap and Tc should be significantly higher and 
we expect that without fluctuations accounted it can be 
around Tc ~ 1.3K. Similarly for optimally doped 100% 
semiconducting SWCNTs the should increase from ob- 
served in [22-23] ~ 15 K to the unsuppressed (by the 
inverse proximity effect of low Tc metallic tubes) Tc of 
19-20 K. The effect of Tc suppression similar to discussed 
here has been observed in alkali metal fulelride molec- 
ular alloys of Ax{CeQ)x{C7o)i-x (24) and adding non- 
superconducting component, i.e. C70 molecules, which 
do not show any superconducting pairing due to symme- 
try reasons and probably due to weaker electron-phonon 
coupling, strongly suppressed Tc from 19 K in 100% Cgo, 
i.e. in KsCqo to Tc=10 K in 20% substituted C70 al- 
loy. The experiments with selectively separated metallic 
and semiconducting SWCNTs, which now become avail- 
able by new methods of effective separation will allow to 
check the validity of presented here simple model and to 
clarify the role of quantum fluctuations, which has not 
been accounted here. 

We introduced a microscopic model of superconduc- 
tivity in a bundle of a mixture of carbon nanotubes. 



We have studied the dependence of a spatially averaged 
superconducting gap < A > on the fraction of semi- 
conducting SWCNT (having higher pairing strength) in 
the bundle at different temperatures. Note that for in- 
homogeneous nanoscale systems the dependence Tc{< 
A(T — 0) >) for different concentration a may be non- 
linear, as a manifestation of the breakdown of the BCS 
theory for bulk materials. Indeed, our calculations of 
Tc{< A(T = 0) >) show a kink at a = 0.5. The reason 
is that the bundle is a highly inhomogeneous system. At 
a < 0.5, below the percolation threshold for a 2D tri- 
angular lattice, the bundle can be seen as a collection 
of finite islands of "good" superconductors (doped semi- 
conducting nanotubes), diluted by normally conducting 
material (metallic nanotubes). Such islands demonstrate 
significantly suppressed superconductivity, even in the 
mean field description, due to the enhanced inverse prox- 
imity effect. Note, that our mean- field BdG model is 
unable to predict and properly describe quantum phase 
fluctuations of the order parameter in quasi one dimen- 
sional systems, where the superconductivity will be sup- 
pressed even stronger. Future research using, for exam- 
ple, Ginzburg-Landau inhomogeneous equations^ i^^i^'^i^^ 
is necessary to describe such kind of effects. Because the 
dynamics of Cooper pairs in doped carbon nanotubes can 
be more close to the diffusive regime, the Usadel equa- 
tions can be applied to calculate the finite conductivity 
at T « Tc^^. 
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